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Abstract 

Background: Caligid copepods, also called sea lice, are fish ectoparasites, some species of which cause significant 
problems in the mariculture of salmon, where the annual cost of infection is in excess of €300 million globally. At 
present, caligid control on farms is mainly achieved using medicinal treatments. However, the continued use of a 
restricted number of medicine actives potentially favours the development of drug resistance. Here, we report 
transcriptional changes in a laboratory strain of the caligid Lepeophtheirus salmonis (Kr0yer, 1837) that is moderately 
(~7-fold) resistant to the avermectin compound emamectin benzoate (EMB), a component of the anti-salmon louse 
agent SLICE® (Merck Animal Health). 

Results: Suppression subtractive hybridisation (SSH) was used to enrich transcripts differentially expressed between 
EMB-resistant (PT) and drug-susceptible (S) laboratory strains of L. salmonis. SSH libraries were subjected to 454 
sequencing. Further L. salmonis transcript sequences were available as expressed sequence tags (EST) from 
GenBank Contiguous sequences were generated from both SSH and EST sequences and annotated. Transcriptional 
responses in PT and S salmon lice were investigated using custom 15 K oligonucleotide microarrays designed using 
the above sequence resources. In the absence of EMB exposure, 359 targets differed in transcript abundance 
between the two strains, these genes being enriched for functions such as calcium ion binding, chitin metabolism 
and muscle structure, y-aminobutyric acid (GABA)-gated chloride channel (GABA-CI) and neuronal acetylcholine 
receptor (nAChR) subunits showed significantly lower transcript levels in PT lice compared to S lice. Using RT-qPCR, 
the decrease in mRNA levels was estimated at ~1. 4-fold for GABA-CI and ~2.8-fold for nAChR. Salmon lice from the 
PT strain showed few transcriptional responses following acute exposure (1 or 3 h) to 200 \ig of EMB, a drug 
concentration tolerated by PT lice, but toxic for S lice. 

Conclusions: Avermectins are believed to exert their toxicity to invertebrates through interaction with glutamate- 
gated and GABA-gated chloride channels. Further potential drug targets include other Cys-loop ion channels such 
as nAChR. The present study demonstrates decreased transcript abundances of GABA-CI and nAChR subunits in 
EMB-resistant salmon lice, suggesting their involvement in avermectin toxicity in caligids. 
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Background 

The development of synthetic insecticides and parasiti- 
cides has revolutionised the control of arthropod pests 
and parasites of agricultural, medical and veterinary 
importance. However, the overuse of chemical control 
agents tends to promote the development of heritable 
insecticide or drug resistance [1], with life-cycle traits of 
targeted organisms often favouring the ability to adapt 
genetically, such as high reproductive potential and short 
generation span. Drug resistance can follow from pre- 
existing mutations in resistance genes, termed resistance 
alleles, which under normal conditions are rare and have 
a negative or neutral effect on fitness but which, under 
conditions of exposure to control agents, afford fitness 
benefits and can become enriched in the target species' 
gene pool given persisting selection pressure from the 
control agent [2]. By the time resistance becomes appar- 
ent as treatment failure, resistance alleles have usually 
already reached high frequencies in the gene pool [2] . 

Caligid copepods, also called sea lice, are common ecto- 
parasites of marine fish [3]. One species, the salmon louse 
(Lepeophtheirus salmonis (Kroyer, 1837)), has emerged as 
a serious problem in mariculture of Atlantic salmon 
{Salmo salar Linnaeus, 1758) in the Northern hemisphere. 
The annual cost of sea louse infection to the global salmon 
farming industry has been estimated at €300 million, with 
the majority of this accounted for through expenses ac- 
crued from treatments with veterinary medicines [4] . Only 
a limited range of anti-sea louse drugs are available and 
licensed for the treatment of fish, and the continued use 
of a relatively small number of compounds creates a situ- 
ation potentially favouring the development of drug resist- 
ance [5]. In the salmon louse, losses of efficacy have been 
reported for a number of control agents including organo- 
phosphates [6], pyrethroids [7], hydrogen peroxide [8] and 
avermectins (AVMs) [9,10]. 

The commonly used anti-sea louse treatment SLICE* 
(Merck Animal Health) contains the avermectin com- 
pound emamectin benzoate (EMB) (Stone et al, 1999). 
SLICE" is administered orally and a one-week treatment 
provides prolonged protection against all host-attached life 
stages of sea lice [11]. Avermectins are also used against 
external and internal parasites of humans and livestock, 
including parasitic nematodes causing the human diseases 
onchocerciasis (River blindness) and lymphatic fUariasis, 
as well as gastrointestinal parasites of sheep, cows and 
horses [12]. The selective toxicity of avermectins against 
ecdysozoan invertebrates is believed to be based on the 
binding and blockage of glutamate-gated (GluCl) and 
y-aminobutyric acid (GABA)-gated (GABA-Cl) chloride 
channels in the invertebrate nervous system [13]. Several 
molecular mechanisms have been suggested as contribut- 
ing factors to the resistance of parasitic nematodes to 
the AVM compound ivermectin (IVM) [14]. Functional 



studies revealed that resistant nematodes can have single 
amino acid mutations in subunits of GluCl and GABA-Cl 
that decrease the channels' sensitivities to the drug [15,16]. 
Furthermore, resistant nematodes may show increased 
expression of ABC (ATP-binding cassette) transporters, a 
group of membrane proteins with members capable of 
mediating the cellular efflux of drugs [17,18]. Finally, aver- 
mectin resistance in insects has been connected to alter- 
ations in drug metabolism [19]. 

Previous studies on potential molecular mechanisms 
of EMB resistance in salmon lice have used the candi- 
date gene approach, i.e. the study of genes that have pre- 
viously been linked to drug resistance in other 
organisms. In particular, such studies have investigated 
salmon louse ABC transporters [9,20] and GABA-Cl and 
GluCl subunits [21]. However, potential outcomes 
resulting from such candidate approaches are limited by 
the current knowledge of the biological process being 
studied. To avoid this limitation, a transcriptomic 
approach was followed in the present study, in which 
microarray analysis was used to compare mRNA re- 
sponses between drug susceptible and moderately 
(~7-fold) EMB resistant laboratory strains of salmon lice. 

Results 

Experimental design 

Two laboratory-maintained strains of salmon lice were 
used in the present study. Strain S is susceptible to all 
currently licensed anti-sea louse treatments including 
EMB, while strain FT is moderately resistant to EMB 
(Figure 1) [9]. Adult male salmon lice were used for the 
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Figure 1 Susceptibility of salmon louse laboratory strains to 
emamectin benzoate (EMB). Toxicity responses of adult male 
salmon lice (Lepeophtheirus salmonis) laboratory strains S and PT in 
24 hours immotility bioassays. Symbols represent the immotility 
response observed in one of duplicate beakers often individuals 
included for each combination of strain and treatment Dose-response 
relationships (solid lines) were established using probit analysis, and 
yielded EC50 values of 73.9 ng (95% confidence intervals: 
58.9 - 92.0 \iq L"^) for the S strain and 6423 (jg L'' (6423 - 9574 |jg L^) 
for the PT strain. 
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transcriptomic analyses as they are considered to pro- 
vide a more steady physiological state than adult 
females, which undergo considerable morphological 
change following fertilisation and which are subject to 
repeated cycles of egg production. Salmon lice were 
collected from host fish anaesthetised in 100 mg L'^ 
2-phenoxyethanol and were then allowed to recover in 
aerated seawater for 2 hours before use. 

In experiment 1, RNA expression profiles were 
analysed in lice sampled directly after recovery in order 
to reveal differences in constitutive gene expression 
between strains. Experiment 2 investigated the effects of 
short-term (1 to 3 hours) exposure to 200 |ig L'^ 
of EMB, compared to seawater and carrier controls 
(Figure 2). This concentration of EMB resulted in >95% 
immotility of S lice after 24 hours, but had no observ- 
able effects on PT lice (Figure 1). We hypothesised that, 
should the tolerances of PT lice to EMB require tran- 
scriptional regulation, this should become apparent in 
transcriptomic profiles during the early stages (1 and 
3 hours) of exposure. In addition to being investigated 
for effects of EMB, the data from experiment 2 were also 
analysed with respect to differences in constitutive tran- 
script expression between strains. For this purpose, data 
for the one and three hour sea water controls from 
experiment 2 were pooled for each strain in order to ob- 
tain the same level of replication as used in experiment 



Experiment 1 



1 (n = 6 biological replicates, i.e. pools of four L. 
salmonis). Control data from experiment 2 were pooled 
between 1 and 3 hour time points as we had previously 
observed that exposure of salmon lice to seawater for up to 
12 hours after collection from the host has no effect on re- 
sponses to EMB in water-borne bioassays (data not shown). 

Custom microarray design 

For transcriptomic analyses, custom Agilent 15 K feature 
60mer oligonucleotide microarrays were designed using 
sequences derived from salmon louse suppression 
subtractive hybridisation (SSH) libraries created in the 
present study together with expressed sequence tags 
(EST) publicly available in GenBank. Although the 
microarrays differed in a minor way between experi- 
ments 1 and 2 due to continuous elaboration/modifica- 
tion, the features included on the microarrays employed 
largely overlap (Table 1). To obtain a pool of salmon 
louse cDNA sequences enriched for transcripts differen- 
tially expressed between the EMB-resistant (PT) and 
drug-susceptible (S) salmon louse strains, two SSH librar- 
ies were constructed, corresponding to subtractions 
between strains in both directions. A pool of both libraries 
was subjected to Roche 454 sequencing, producing a total 
of 94,834 reads (N50 value of 289 nucleotides). The 
assembly of contiguous sequences (contigs) from sequence 
reads provided 1,916 annotated (BLASTx e-value <10'*) 



Strain PT 
(Resistant) 



untreated 



n = 6 



Strain PT 
(Resistant) 



Strain S 
(Susceptible) 



untreated 



n = 6 



Experiment 2 



Strain S 
(Susceptible) 



Seawater 1 h PEGg^^l h 



Seawater 3 h 



n = 3 



EMB 1 h 



Seawater 1 h PEGg^^l h 



EMB3h 



Seawater 3 h 



n = 3 



EMB 1 h 



EMB3 h 



Figure 2 Design of microarray experiments. Microarray analysis was used to measure transcriptomic profiles in drug-susceptible (S) or EMB 
resistant (PT) salmon louse strains. In experiment 1, adult male salmon lice were collected for RNA extraction directly after removal from host fish. 
In experiment 2, salmon lice were removed from hosts and then subjected to exposure of seawater containing 200 (jg L'' of EMB for 1 or 
3 hours before sampling. Control treatments included exposures to seawater, or seawater containing the solvent PEG300, which was used to 
solubilise EMB. The number of pools of four salmon lice sampled for each combination of treatment and strain is given as n. 
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Table 1 Composition of features on custom L. salmonis 



oligo microarrays 



Probe type 


Experiment 1 
(AMADID# 039612) 


Experiment 2 
(AMADID # 033382) 


Annotated 


1 3, 542 


10,055 


Unannotated 


1 ,566* 


5,052* 


Control probes 


100 


100 


Agilent controls 


536 


536 


Total 


15744 


15744 



* Oligo probes were designed to bothi the forward and reverse complement 
sequence for all targets that were not suitably annotated but included on the 
custom microarrays. 



and 783 un-annotated target sequences. KEGG functional 
analysis of the annotated sequences revealed a large repre- 
sentation (53%) of genes involved in metabolism (data not 
shown). In addition, 129,225 L. salmonis ESTs (> 100 bp) 
were obtained from GenBank and assembled into contigs, 
providing a further 10,056 annotated (BLASTx e-value 
<10'*) and 2,526 un-annotated target sequences for the 
design of oligo probes to be used in the microarray designs 
(Table 1). 

Analysis of strain differences in constitutive gene expression 

To determine constitutive differences in gene expression 
between the PT and S strains, mRNA expression profiles 
were analysed in adult male salmon lice sampled in 
February 2012 (experiment 1) and May 2011 (seawater 
controls of experiment 2). When data from each experi- 
ment were analysed including only features present on 
both microarrays, similar numbers of features were 
found to be differentially expressed between strains in 
experiment 1 and experiment 2 (1,113 and 1,280 features 
respectively; Figure 3). Comparison of these two feature 
lists revealed that only 359 features were reported as 
being significantly differentially expressed between 
strains in both experiments. Of these, 294 (82%) showed 
the same direction of strain differences in the two exper- 
iments (Figure 3) and represented 226 genes of which 
57% were annotated. These 226 genes were arranged by 
significance of the expression differences determined in 
experiment 1. 

Genes that were represented in the top 100 most sig- 
nificantly differentially expressed transcripts are detailed 
in Table 2, which includes a substantial number of cyto- 
skeleton proteins (26%) and proteases (12%). Enrichment 
analysis of the 294 features resulting from comparison of 
expression between strains was performed with respect 
to the gene ontology (GO) annotation representation on 
the microarray. Nine GO attributes were found to be 
significantly over-represented (Additional file 1), with 
calcium ion binding, structural constituent of muscle 
and actin binding being shown to be the most signifi- 
cantly over-represented GO terms. 



To confirm findings from microarray analyses, tran- 
script abundance was analysed for a sub-set of signifi- 
cantly differentially expressed genes using RT-qPCR. 
Genes were selected on the basis of potential signifi- 
cance as pharmacological targets of EMB (GABA-Cl 
subunit alpha and neuronal acetylcholine receptor sub- 
unit a3), or detoxification mechanisms (cytochrome 
P450 isoforms, carboxylesterase). Maltase-glucoamylase 
was further included because of its high level of differen- 
tial expression (105-fold) between salmon louse strains. 
RT-qPCR analysis found that transcripts of nAChR a-3 
were ~3. 1-fold and ~2.6-fold less abundant in the PT 
than the S strain in experiments 1 and 2 respectively, 
which confirmed trends observed in the microarray ana- 
lyses (Table 3). Similarly, RT-qPCR demonstrated that 
PT lice showed significantly lower levels of GABA-Cl 
a-subunit mRNA expression compared to the S strain 
(1.4-fold and 1.6-fold in experiments 1 and 2, respect- 
ively; Table 3), although differences were marginal. 
Isoforms of cytochrome P450 and carboxylesterase, i.e. 
enzymes potentially involved in detoxification, were 
found to show higher mRNA expression levels in the PT 
compared to the S strain in experiment 1, but not 
experiment 2. Transcript levels of maltase-glucoamylase 
were much lower in PT than S lice in experiment 1, 
whereas in experiment 2 the mRNA expression was 
moderately increased in PT compared to S lice (Table 3). 

Effects of short-term exposure to EMB on transcript 
profiles in salmon lice 

Three-way ANOVA of the microarray expression data 
from experiment 2 showed that a large proportion (55%) 
of the total number of features that passed quality filter- 
ing (n = 10,804) was affected by the factor strain. In con- 
trast, the factor treatment had a comparatively small 
influence on gene expression (Table 4). To confirm these 
microarray results, transcript abundances were deter- 
mined by RT-qPCR for six selected genes. Genes were 
selected so as to include a number of qualitatively differ- 
ent expression profiles, detectable fold-changes and 
selected candidate genes (GABA-a and AChR-a3 
subunits) (Figure 4). A high degree of correlation was 
observed between expression values measured by both 
methods (Pearson correlation coefficients (r) of 0.71 to 
0.99; p <0.0001). 

To further investigate the effects of EMB exposure, a 
list of those features that were significantly affected by 
EMB treatment or for which significant interactions be- 
tween treatment and other factors were observed (treat- 
ment X strain; treatment x time; treatment x strain x 
time) was compiled. This list comprised a total of 2,020 
features, of which at least 35% were involved in metabol- 
ism (KEGG functional classification, data not shown). 
Transcript abundance profiles for the 2,020 features 
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Experiment 1 




Experiment 2 



Experiment 1 
Expression 
change (PT/S) 



Experiment 2 
Expression 
cliange (PT/S) 



Number of 
features 



Number of 
genes 



Down 
Up 

Down 
_Up 



Down 

Up 

Up 

Down 



235 (66%) 
59(16%) 
40(11%,) 
25 (7%) 



180 n 

46 J 

37 

21 



226* 



Total 



359 



284 



Figure 3 Comparison of microarray features differentially expressed across salmon louse strains between experiments. Genes 
differentially expressed between two salmon lice strains differing in EMB susceptibility (S, PT) were determined by microarray analysis in two 
independent experiments. Features that showed differential expression between strains (fold change > 1.3, p < 0.05) were derived for both 
experiments. The Venn diagram only includes features studied in both experiments. *A total of 226 genes showed comparable differential 
expression between salmon louse strains in two independent microarray experiments and were therefore selected for further analysis of 
biological function. 



were further subjected to network clustering using the 
BioLayout Express""^ application [22]. This resulted in 
the resolution of 59 clusters with a minimum cluster size 
of four features. The two main clusters 1 and 2 
contained 418 and 62 features, respectively, that showed 
fold changes > 1.3 across all conditions (Figure 5A & 
B). Within the two clusters, expression profiles were 
characterised by pronounced responses in S lice follow- 
ing 1 h of EMB exposure (down-regulation in cluster 1, 
up-regulation in cluster 2), and a full or partial return to 
basal expression levels after 3 h of EMB exposure. More- 
over, for genes in both clusters, few if any responses to 
EMB exposure were observed in the PT strain. 

The most significantly differentially regulated genes 
from cluster 1, and the responses to EMB exposure in 
both salmon louse strains are summarised in Additional 
file 2. An analogous selection of genes from cluster 2 is 
provided in Additional file 3. Genes in cluster one in- 
clude the GABA-Cl and nAChR a-3 subunits which are 
potential targets for EMB, with nAChR a-3 also having 
been identified as constitutively differentially expressed 



between salmon louse strains. Genes in cluster two in- 
cluded a glutathione-S -transferase isoform and a nAChR 
a-3 precursor. Enrichment analysis of the features in 
cluster one is detailed in Additional file 4, showing sig- 
nificant over representation of twelve and under repre- 
sentation of four GO attributes. Chitin binding, calcium 
ion binding and hydrolase activity were the most signifi- 
cantly over-represented attributes and nucleic acid bind- 
ing was identified as the most significant under 
represented attribute. 

Discussion 

Using transcriptional profiling in comparative studies of a 
drug-susceptible and an EMB-resistant salmon louse 
strain, this study demonstrated the reduced constitutive 
mRNA expression of subunits of certain ligand-gated ion 
channels (LGIC) in the EMB resistant strain, namely a 
GABA-gated chloride channel subunit (GABA-Cl, -1.4- 
fold decreased) and a neuronal acetylcholine receptor 
subunit (nAChR a-3; ~2.8-fold decreased). The toxicity of 
EMB and other AVMs against ecdysozoan invertebrates is 
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Table 2 Genes showing differential expression between salmon louse strains differing in EMB susceptibility 


Accession no. 


Annotation 


Experiment 1 




Experiment 2 








Fold change (PT/S) 


p-value 


Fold change (PT/S) 


p-value 


Ligand gated ion channel (<1%) 










ADD24187.1 


Neuronal acetylcholine receptor subunit alpha-3 


-3.21 


7.59E-05 


-2.52 


4.29E-03 


Metabolism of xenobiotics (<l%) 










ACOl 5001.1 


Cytochrome P450 3A24 


-1.48 


1 .95E-03 


-1.31 


4.13E-02 


Regulation of synapse development (<l%) 










ADD24462.1 


Cerebellin-3 


-1.86 


3.60E-04 


-4.91 


3.53E-06 


Eye lens proteins (<1 


'%) 










ADD38in.l 


Beta-crystallin Al 


2.33 


6.79E-05 


1.37 


1.35E-02 


Transporters (2%) 












NP_001 116712.1 


Solute carrier family 8 (sodium/calcium exchanger) 


-2.10 


3.63E-04 


-4.38 


1 .74E-05 


EFX88361.1 


Alpha subunit of Na+/K+ ATPase 


-1.67 


2.21 E-04 


-2.08 


5.34E-05 


ACOl 261 3.1 


Excitatory amino acid transporter 3 


-1.42 


1.10E-04 


-1.34 


4.96E-02 


Cytoskeleton proteins (26%) 










ADD38332.1 


Troponin 1 


-1.77 


4.85 E-05 


-2.61 


1.17E-06 


XP_001 950563.1 


Muscle LIIVl protein 


-1.67 


1 .57E-04 


-2.30 


1.32E-06 


AAA1 7371.1 


Fast myosin heavy chain 


-2.11 


1 .73E-03 


-6.23 


1 .68E-06 


ACOl 3 186.1 


Myosin light chain alkali 


-2.65 


1 .69E-05 


-3.38 


8.46E-06 


ACOl 2924.1 


Myosin light chain alkali 


-2.06 


1.16E-05 


-3.69 


9.66E-06 


ABU41018.1 


Collagen alpha-1 chain 


-1.38 


1.51 E-03 


-6.85 


2.61 E-05 


EFV6 1840.1 


Smoothelin 


-1.43 


4.03E-04 


-1.85 


5.25 E-05 


ACOl 2887.1 


Troponin T 


-1.73 


1 .04E-03 


-2.93 


8.12E-05 


ADV40202.1 


Troponin 1 


-1.94 


1 .87E-03 


-2.93 


3.23E-04 


AAW22542.1 


Myosin light chain 


-1.76 


1 .OOE-03 


-2.23 


3.73E-04 


ACOl 0528.1 


Troponin C, isoform 1 


-1.91 


1 .69E-03 


-2.22 


5.12E-04 


ACOl 2630.1 


Troponin C, isoform 1 


-1.61 


3.01 E-05 


-1.51 


5.92E-04 


ACOl 2421.1 


Tropomodulin 


-1.51 


1.17E-04 


-2.37 


1 .01 E-03 


ACOl 2794.1 


Troponin C, isoform 1 


-2.07 


6.14E-05 


-1.60 


2.85E-03 


ACOl 475 1.1 


Troponin C, isoform 1 


-2.12 


4.94E-05 


-1.56 


4.18E-03 


ACOl 1818.1 


Torso-like protein precursor 


-1.83 


5.11 E-05 


-1.82 


9.67E-03 


ACOl 1077.1 


Troponin C, isoform 1 


-1,92 


6.27E-05 


-1.48 


1 .56E-02 


Regulation of actin cytoskeleton (<l%) 










XP_002407362.1 


Paxillin 


-2.79 


8.32E-07 


-1.48 


2.28E-03 


Calcium transport (<l%) 










NP_001 03271 9.1 


Sarco/endoplasmic reticulum calcium transporting ATPase 


-2.15 


8.58E-07 


-2.46 


Z15E-04 


Calcium binding (1%) 










XP_002734090.1 


Calmodulin-like 


-1.60 


641 E-05 


-2.46 


3.25E-04 


ACOl 1757.1 


Sarcoplasmic calcium-binding protein, beta chain 


-3.18 


2.05 E-06 


-3.58 


1.28E-03 


Cuticle proteins (4%) 












ADD24515.1 


Cuticle protein 6 


-2.29 


1 .79E-03 


-12.28 


4.75E-04 


ABU41 025.1 


Cuticle protein 


-2.57 


2.78E-04 


-12.91 


1.88E-03 


ACOl 4885.1 


Cuticle protein CP14.6 precursor 


-1.77 


7.50E-04 


-3.40 


4.67E-03 
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Table 2 Genes showing differential expression between salmon louse strains differing in EMB susceptibility (Continued) 



Proteolysis (12%) 












ADD38666.1 


Matrix metalloproteinase-9 


-2.29 


1 .07E-05 


-7.35 


1 .84E-06 


ADD38283.1 


Kunitz/BPTI-lil<e toxin 


-2.45 


543E-06 


-1.79 


2.50E-04 


ABU41 053.1 


Metalloproteinase 


-5.61 


7.00E-06 


-4.38 


3.99E-04 


BAG74353.1 


Metalloproteinase 


14.83 


8.20E-05 


18.61 


4.64E-04 


AC01 1096.1 


Serine carboxypeptidase CPVL precursor 


1.58 


2.44E-04 


1.36 


1.21 E-02 


ABU4in7.1 


IVletalloproteinase 


-2.72 


6.52E-05 


-1.41 


1.82E-02 


AAS91 793.1 


Intestinal trypsin 2 precursor 


7.73 


1 .90E-06 


1.82 


2.52E-02 


AAS91 795.1 


Intestinal trypsin 4 precursor 


1.48 


6.69E-04 


1.42 


2.79E-02 



Annotated genes (57%) were sorted by the significance of differential expression between strains in experiment 1 and arranged by biological function. Features 
withi identical annotation were removed prior to categorising biological function. Also indicated is the percentage representation for each functional category in 
the total number of annotated genes (129 genes). 



reported to be based mainly on their interaction with 
another class of LGICs, the glutamate-gated chloride 
channels (GluCl) [23-25] although GABA-Cls are also be- 
lieved to be pharmacological targets of AVMs [26] . While 
nAChRs are traditionally not considered to be implicated 
in the toxic action of AVMs in ecdysozoans, they can be 
allosterically modulated by IVM [27]. This study's finding 
that mRNA levels of GABA-Cl and nAChR subunits are 
decreased in EMB-resistant salmon lice suggests that these 
LGICs may represent potentially additional target sites for 
AVMs in sea lice. When transcriptomic profiles of EMB- 
resistant and -susceptible laboratory salmon louse strains 
were investigated in response to short term (1-3 h) aque- 
ous EMB exposure, a number of transcriptional responses 
to the treatment were observed in the drug- susceptible S 
lice, but few responses were found in the EMB-resistant 
PT strain. While we cannot exclude the possibility that 
EMB exposure might provoke more pronounced 
transcriptomic responses in PT lice at later time points, 
differential toxicity of EMB between the louse strains is 



apparent as early as 5 hours post-exposure (data not 
shown). This suggests that the mechanism of resistance, 
whether constitutive, responsive or a combination of both, 
must have occurred before this time point. The present 
study only considered levels of mRNA. Molecular mecha- 
nisms underlying differential susceptibilities between the 
studied strains could also include post-transcriptional 
regulatory mechanisms, such as mRNA processing and 
degradation, translation and protein degradation [28], 
however, these lay outside the scope of the present study. 

AVMs are used against parasitic and pest species of 
ecdysozoan invertebrates including nematodes, insects 
and mites, and more recently against crustaceans. While 
the molecular target sites of AVMs in crustaceans are 
unknown, GluCls are generally considered to be the 
main pharmacological targets of IVMs in nematodes and 
insects [29,30]. The GluCls form an invertebrate-specific 
subgroup of the large 'Cys-loop' subfamily of LGIC. 
'Cys-loop' LGICs have a pentameric structure, and are 
composed of either the same type of subunits or two to 



Table 3 Gene expression measured by RT-qPCR in salmon lice from two strains differing in EMB susceptibility 

Accession no. Annotation Experiment 1 Experiment 2 



Microarray RT-qPCR Pearson Microarray RT-qPCR Pearson 

correlation correlation 







p-value 


Fold 
change 
(PT/S) 


p-value 


Fold 
change 
(PT/S) 


r 


p-value 


p-value Fold 
change 
(PT/S) 


p-value 


Fold 
change 
(PT/S) 


r 


p-value 


ADD24187.1 


Neuronal acetylcholine 
receptor subunit a3 


7.59E-05 


-3.21 


0.000 


-3.08 


0.93 


<0.0001 


429E-03 


-2.52 


0.009 


-2.64 


0.89 


0.0001 


EFN73916.1 


GABA receptor subunit 
alpha 


NS 


-1.19 


0.015 


-1.35 


N/A 


N/A 


2.05 E-04 


-1.75 


0.005 


-1.56 


0.89 


0.0001 


XP_003494528.1 


Cytochrome p450 18al 


3.42E-08 


3.34 


0.000 


3.14 


0.95 


<0.0001 


N/A 




NS 


1.00 


N/A 


N/A 


AAS 13464.1 


Cytochrome p450 ISal 


8.93E-04 


2.14 


0.000 


1.82 


0.95 


<0.0001 


N/A 




NS 


1.09 


N/A 


N/A 


NP_001 1361041 


Carboxylesterase 


1.13E-04 


1.38 


0.041 


121 


0.65 


0.02 


NS 


1.12 


NS 


1.12 


N/A 


N/A 


XP_797271.2 


Maltase-glucoamylase 


9.12E-05 


-104.50 


0.000 


-19485 


0.99 


<0.0001 


7.10E-03 


1.72 


0.017 


2.00 







Abbreviations: NS Not significantly different (p<0.05), N/A Not applicable. 

For selected genes, mRNA expression was measured in salmon lice from two strains differing in EMB susceptibility (see Figure 1 for information on strains), results 
compared to findings obtained in microarray analyses. Significance (p < 0.05) assessed by f-test (Welch) for microarray analysis and one-way ANOVA for RT-qPCR 
analysis. Fold changes > 1.3 are underlined. 
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Table 4 Number of differentially expressed features 
identified in microarray analysis of samples from 
experiment 2 


Experimental factor 


Number of features 


Strain 


5940 (55%) 


Treatment 


369 (3%) 


Time 


968 (9%) 


Strain x Treatment 


406 (4%) 


Strain x Time 


950 (9%) 


Time X Treatment 


1309 (12%)) 


Strain x Treatment x Time 


1701 (16%)) 



Significance (p-value < 0.05) assessed by three-way ANOVA. Also indicated is 
the percentage of the total number of features (10, 804) used in the analysis. 



three different subunit isoforms. Consistent with the role 
of GluCl as the main target of AVMs, IVM-resistant 
strains of invertebrates can show mutations changing 
the expression levels or the peptide sequence of channel 
subunits [23,24]. A GluCla subunit has been cloned in 
L. salmonis [21] and while GluCla was represented 
amongst the microarray targets used for this study, no 
difference in mRNA expression was observed between 
salmon lice of the two studied strains, or amongst those 
subjected to control and sublethal EMB treatments (data 
not shown). 

Apart from GluCl, further LGICs are known to inter- 
act with AVMs. For instance, IVM modulates the activity 
of nematode GABA-Cl [26], and can exert directly acti- 
vating or potentiating effects on vertebrate glycine-gated 
chloride channels [31]. Moreover, AVMs can modulate 
the activity of cation-LGICs such as the a-7 nAChR [27] 
and the ATP-gated P2X4 receptors [32]. A number of 
observations involving drug-resistant insects and nema- 
todes support the hypothesis that LGICs other than 
GluCl constitute further toxicologically relevant targets 
of AVMs in invertebrates. Cyclodiene-resistant fruit flies 
having a single amino acid mutation in a GABA-Cl 
showed a moderate degree of cross-resistance to IVM 
[24]. A null mutation in a histamine-gated chloride 
channel also conferred moderate IVM resistance in 
Drosophila melanogaster Meigen, 1830 [33], and a novel 
dopamine-gated ion channel (HcGGR3) was significantly 
down-regulated in an AVM-selected strain of the nema- 
tode Haemonchus contortus (Rudolphi, 1803) [34]. The 
observation in this study that EMB-resistant salmon lice 
show decreased mRNA levels of nAChR and GABA-Cl is 
consistent with findings in the literature cited above, and 
suggests a role for nAChR and GABA-Cl as additional 
pharmacological targets of EMB in salmon lice. It is 
worth noting in this context that observed changes in 
nAChR expression could also relate to previous exposure 
of PT lice to compounds interfering with cholinergic 



neurotransmission such as the organophosphate (OP) 
anti-sea louse drug azamethiphos [35]. However, decreases 
in nAChR expression are not among typical molecular re- 
sponses associated with OP resistance in insects [1,36]. 
While decreased expression of nAChR has been observed 
in nematodes resistant to imidazothiazoles and other drugs 
that directly target nAChR channels [14], these are classes 
of compound that have not been used against sea lice. 

Apart from modifications of the molecular targets, 
biocide resistance in pests and parasites can result from 
increased elimination of the chemical as a result of 
up-regulation of biotransformation enzymes and/or drug 
transporters. The superfamily of cytochrome P450s 
(CYPs) contains heme-thiolate proteins that function as 
monooxygenases, many of which are involved in drug 
metabolism [37,38]. CYPs play important roles as chem- 
ical resistance genes in insects [1,39], but their roles in 
the biochemical defence against toxicants in crustaceans 
are less well understood. The microarrays used in the 
present study included probes representing a number of 
target sequences that are annotated as CYPs (experiment 
1: 18 probes; experiment 2: 14 probes). These partial 
L. salmonis CYP sequences could not be unequivocally 
attributed to specific CYP families, as this would require 
establishing the overall degree of amino acid similarity 
to CYP family members [40]. However, based on 
BLASTx annotation (e-values of < 10' ), most of these 
sequences could be provisionally allocated to CYP clans, 
which are higher-order groupings that combine phylo- 
genetically related CYP families [40] (Clan 2: 6 
sequences; clan 3: 5 sequences; mitochondrial clan: 3 
sequences). Two targets showing similarity to CYP15A1 
and CYP18A1 (both clan 2) differed in mRNA expres- 
sion in L. salmonis from experiment 1, but not experi- 
ment 2. Moreover, compared to the EMB susceptible 
strain a target showing similarity to CYP3A24 (clan 3) 
had lower expression levels in the EMB resistant strain, 
with a moderate (1.48-fold) difference in transcript 
abundance found between the strains. In insects, 
CYP15A1 is involved in juvenile hormone synthesis [41], 
whereas CYP18A1 functions to inactivate ecdysteroids 
[42]. Crustaceans possess homologues to both these clan 2 
CYPs [43], but littie is known of their function. In the 
green shore crab {Carcinus maenas) expression levels of 
two CYPs from clan 2 were affected by both the moulting 
cycle and previous exposure to xenobiotics [44]. The 
differences in CYP mRNA levels found in this study were 
relatively small and/or variable between experiments and 
therefore do not provide clear evidence for an involvement 
of CYPs in the differential EMB susceptibility found in the 
salmon louse strains that were studied. 

Carboxylesterases are another class of enzymes that 
can confer insecticide resistance [1]. In this study, the 
expression of one carboxylesterase was moderately 
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Figure 4 Gene expression responses at early time points of exposure to EMB. Shown are the relative expression ratios (RER) measured by 
microarray and RT-qPCR analysis of S and PT salmon louse strains after 1 and 3 hours exposure to 200 [ig l'^ EMB, seawater (SW), or the solvent 
PEG300 (SOL) . Data are Log2 RER ± SE (n = 3). 
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enhanced in EMB-resistant salmon lice in experiment 1, 
but no significant differences in expression were ob- 
served between strains in experiment 2. Accordingly, the 
data provides no evidence for a role of carboxylesterase 
in EMB resistance of salmon lice. 

ABC (ATP-binding cassette) proteins are a family of 
membrane-bound transporters mediating the transport 
of a diverse array of substrates across biological mem- 
branes [45]. Certain ABC proteins are drug efflux trans- 
porters located in the cell membrane, and have roles in 
the biochemical defence against toxicants [46]. The ABC 



transporter P-glycoprotein transports IVM [47] and has 
relevance as a biochemical factor limiting the drug's 
toxicity in mice and nematodes [18,48]. It has been 
suggested that P-glycoprotein could be implicated in the 
resistance of pests and parasites to AVMs [49-51], and a 
salmon louse homologue of P-glycoprotein called SL- 
PGYl has been cloned [9]. In the present study, mRNA 
expression of SL-PGYl was unaffected by EMB treat- 
ment and did not differ between S and PT strains (data 
not shown), confirming similar data previously reported 
for these strains [9]. 
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In addition to comparing transcriptional profiles be- 
tween louse strains in the absence of drug exposure, the 
present study also addresses transcriptomic responses of 
both louse strains to short-term exposures (1 to 3 hours) 
to 200 |ag L"' of EMB. The selected EMB concentration 
is not toxic to the resistant salmon louse strain, but in 
the drug-susceptible strain it results in close to 100% 
immotility after 24 hours exposure. After 1 to 3 hours of 
exposure, few transcriptomic responses were observed in 
the resistant strain, and therefore these data do not 
provide evidence for the presence of specific resistance 
mechanisms requiring transcriptional regulation. In the 
drug-susceptible strain, a complex array of transcriptional 
responses to EMB was observed. This finding was 
expected, as this level of EMB exposure was known to 
result in significant toxic responses after 24 hours. Early 
transcriptional responses are therefore likely to comprise 
both adaptive and general stress responses, as well as tran- 
scriptional changes reflecting the toxic action of the drug. 

Conclusions 

The transcriptional profiling of a drug susceptible and 
an EMB resistant laboratory strain of salmon lice re- 
vealed reduced mRNA expression of a number of LGIC 
subunits in the EMB resistant lice. GluCl, considered to 
be the major target site of AVM drugs in invertebrates, 
displayed similar mRNA expression levels in EMB- 
resistant and reference strains. In contrast, subunits of 
GABA-Cl and nAChR showed decreased mRNA abun- 
dances in the EMB resistant compared to the reference 
strain. While GABA-Cl is considered a secondary target 
of AVMs in invertebrates, nAChR is not traditionally 
considered a target site for AVMs, even though it has 
been shown to interact with AVMs in vertebrates. It is 
possible that nAChR and GABA-Cl represent additional 
EMB target sites in salmon lice, and that the down- 
regulation of these channel subunits in this EMB- 
resistant strain could be related to the resistance pheno- 
type. In the present study no changes were seen in the 
expression levels of biotransformation enzymes and drug 
transporters, both of which classes have been suggested 
to contribute to AVM resistance in other species. Fur- 
ther studies are needed to investigate potential relation- 
ships between the transcriptional changes observed and 
the susceptibility phenotype. 

Methods 

Salmon louse strains 

Two laboratory-maintained salmon louse (L. salmonis) 
strains differing in susceptibility to EMB [9], were used 
in this study. The susceptible strain S was established 
in 2003 using salmon lice from a Scottish farm site 
where no chemical control agents other than hydrogen 
peroxide had been used. The moderately EMB-resistant 



salmon louse strain PT was established in December 
2008 using salmon lice from another Scottish production 
site with reports of variable EMB treatment efficacies. 
The strains have since been cultured under identical la- 
boratory conditions, as described in detail elsewhere [9]. 
In brief, salmon lice were maintained on Atlantic salmon 
[S. salar) with an initial weight of 500-1000 g in circu- 
lar tanks supplied with fresh seawater at ambient 
temperature, using a photoperiod corresponding to nat- 
ural day length. To propagate salmon louse cultures, egg 
strings were allowed to hatch and develop to copepodids, 
which were used to inoculate a tank containing fresh host 
fish. Prior to the collection of salmon lice from hosts, fish 
were anaesthetised with 100 mg L'^ 2-phenoxyethanol. In- 
fection rates were maintained at levels that were unlikely 
to compromise fish welfare. All laboratory infections were 
carried out under UK Home Office licence and appropri- 
ate veterinary supervision. 

Salmon louse exposure experiments 

Adult male salmon lice were collected from anaes- 
thetised host fish as described above and allowed to re- 
cover for 2 hours in aerated filtered seawater at ambient 
sea temperature. To analyse transcript expression in 
salmon louse strains S and PT in the absence of drug 
exposure (microarray experiment 1), adult males were 
collected and preserved in an RNA stabilisation solution 
(4.54 M ammonium sulphate, 25 mM trisodium citrate, 
20 mM EDTA, pH 5.4) prior to storage at -70°C. For 
both strains six pooled samples comprising four individ- 
uals each were taken. The effects of acute seawater- 
borne exposure to EMB on transcript expression were 
investigated in microarray experiment 2. In order to 
establish the EMB concentration for this experiment, 
24 hour bioassays [52] were carried out on salmon louse 
strains S and PT. EMB (technical grade, a gift from 
Merck Animal Health) was solubilised in seawater with 
PEG300 (final concentration 0.01% (v/v)). EMB concen- 
trations used in bioassays were 32.5, 75, 150, 300 and 
600 \xg L'^ with S lice, and 200, 400, 800, 1,600 and 
3,200 |ig L'^ with PT lice. Duplicate glass dishes 
containing 10 salmon lice and 200 mL of exposure solu- 
tion were used per EMB concentration, control (sea- 
water) or solvent control (seawater with 0.01% (v/v) 
PEGsoo)- At the end of 24 hours of exposure, salmon lice 
were recorded as normally motile or immotile upon 
visual examination and stimulation with a fine brush. To 
generate samples for microarray experiment 2, salmon 
lice of either strain were subjected to short term (1 and 
3 hour) exposures to 200 i^g L'^ EMB, a concentration 
that would result in >95% immotility in S lice after 
24 hours but have no effects in PT lice. In addition, the 
experiment comprised seawater and solvent (0.01% (v/v) 
PEG300) controls. For each combination of strain. 
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exposure period and treatment, three pooled samples 
consisting of four salmon lice each were collected for 
later RNA extraction as above. None of the treatments 
had effects on louse motility. At the end of the experi- 
ment (3 hours), water samples were taken and sent to a 
commercial laboratory (Eclipse Scientific, Chatteris, UK) 
for EMB residue analysis (liquid chromatography with 
detection by MSMS). The measured EMB concentration 
in the nominal 200 i^g L'^ EMB treatment was 99.5 ± 
5.2 |ig L'^ EMB. This depletion of solubilised active 
ingredient may be attributed to EMB adsorption to the 
glass containers used for exposure assays [53]. 

RNA extraction and purification 

In microarray and RT-qPCR experiments, samples were 
pools of four adult male salmon lice. Frozen samples 
were ground in liquid nitrogen using a pestle and mor- 
tar, and total RNA was immediately extracted from 
the homogenised sample using TRI Reagent* (Sigma- 
Aldrich, UK), following the manufacturer's protocol. 
After phase separation, RNA was precipitated from the 
aqueous phase by addition of 0.25 volumes isopropanol 
and 0.25 volumes of a high salt buffer (0.8 M trisodium 
citrate; 1.2 M sodium chloride), as recommended for 
samples with high polysaccharide content [54]. The total 
RNA was resuspended in nuclease-free water and fur- 
ther purified using RNeasy columns (Qiagen, UK). For 
the construction of subtracted cDNA libraries, total 
RNA from 60 untreated adult males from either strain 
(S or PT) were pooled and subjected to poly (A)+ RNA 
isolation using the Poly (A) Purist™ kit (Ambion*, UK). 
UV spectroscopy (NanoDrop ND-1000, Thermo Scien- 
tific, USA) was used to confirm purity of the RNA sam- 
ples and establish concentrations, whereas RNA 
integrity was assessed by agarose gel electrophoresis and 
ethidium bromide staining. 

Subtracted cDNA library construction and sequencing 

Suppression subtractive hybridisation (SSH) was used to 
prepare cDNA libraries enriched in transcripts differen- 
tially expressed between strains S and PT using com- 
mercial methods (PCR-Select™, Clontech, Takara Bio 
Inc., USA) . Subtractions were performed in both direc- 
tions, i.e. using cDNA derived from each strain (S or 
PT) either as the tester or the driver. A pool of cDNA 
from each subtraction, containing an equal amount of 
both subtracted cDNA libraries, was used for generating 
a 454 sequencing library using the GS FLX Titanium 
Rapid Library Preparation kit (Roche Applied Science, 
UK), following manufacturer's instructions. Adaptive 
Focus Acoustics'" (AFA™) using the S220 High Perform- 
ance Ultrasonicator (Covaris" Inc., KBiosciences, UK) 
was employed to randomly shear the cDNA, blunt ends 
were repaired and MID adapters ligated to the DNA 



fragments prior to sequencing using the Genome 
Sequencer™ (GS) Titanium FLX instrument (Roche Ap- 
plied Science, UK) (EBI Sequence Read Archive (SRA) 
study ERP002190). GS FLX Titanium library preparation 
and sequencing was performed by The GenePool Gen- 
omics Facility (University of Edinburgh, UK). Sequence 
reads were assembled using the GS De Novo 
Assembler (Newbler) v2.5.3 software (Roche Applied 
Science, UK) using default parameters after trimming of 
MID adapter and primer sequences. 

Sequence assembly from existing L salmonis EST 
resource 

A total of 129,225 sequences (> 100 bp) were downloaded 
during December 2010 for L. salmonis from the GenBank 
EST database (as described by Yasuike et al.) [55], and 
assembled into contigs using default assembly settings of 
The Gene Indices Clustering Tools (TGICL), obtained 
from the Computational Biology and Functional Genom- 
ics Laboratory (The Gene Index Project, USA). Prior to 
sequence assembly, vector sequences were removed using 
SeqMan II 6.1 (DNAStar Inc., USA). 

Salmon louse microarray design 

The assembled contig sequences were annotated using 
BLASTx (Basic Local Alignment Search Tool) searches 
against the non-redundant proteins (nr), UniprotKB/ 
Swiss-Prot (Swissprot) and Reference Proteins (Refseq_ 
Proteins) GenBank databases at the National Centre for 
Biotechnology Information (NCBI), with an annotation 
hit having an expectation value (e-value) of <1 x 10' 
being considered significant. All sequences were fur- 
ther annotated with GO identifiers using Blast2Go soft- 
ware for Windows* using Java Webstart (Centro de 
Investigacion Principe Felipe, Spain). Oligonucleotide 
probes (60mers) were designed to target contig se- 
quences using the eArray Gene Expression (GE) probe 
design tool (Agilent Technologies, UK), employing the 
base composition and best probe methodologies, and 
designed in sense orientation with 3' bias. For each se- 
quence without a significant BLASTx based annotation 
two probes were selected; designed to both forward and 
reverse complement sequences. Standard expression 
microarrays were designed using the eArray custom 
microarray design wizard (Agilent Technologies, UK) for 
an 8 X 15 K design format. Each microarray comprised 
15,744 features including 536 obligatory controls 
(Table 1). An initial design was used for experiment 1 
interrogations (Agilent AMADID No 039612; EBI 
ArrayExpress design A-MEXP-2285). This design incor- 
porated probes designed to target 2,699 sequences that 
were identified when sequencing the subtracted cDNA 
libraries enriched for transcripts differentially expressed 
between the EMB-resistant (PT) and drug susceptible 
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(S) salmon louse strains. Two probes were designed for 
each of the SSH targets. Experiment 2 employed a 
modified design (Agilent AMADID No 033382; EBI 
ArrayExpress design AMEXP-2284). The array designs 
shared 10,251 identical features. 

Microarray analyses 

Labelling protocols are described in detail elsewhere [56]. 
Briefly, for each test sample 250 ng total RNA was used as 
template for the amplification of antisense RNA with the 
incorporation of the modified nucleotide 5-(3-aminoallyl)- 
UTP (aaUTP) into the amplified RNA (aRNA) during the 
in vitro transcription step (TargetAmp™ Aminoallyl-aRNA 
Amplification Kit 101; Epicentre', Cambio Ltd. UK). A 
common reference pool was created through pooling 
equal amounts of all aRNA test samples to be used in the 
experiment. The individual test samples were labelled with 
cyanine 3 (Cy3) and the common reference pool labelled 
with Cy5 mono-reactive dye (GE Healthcare, UK) in dye 
coupling reactions. Unincorporated dye was removed by 
column purification (lUustra Autoseq"" G-50 spin columns; 
GE healthcare, UK), and then dye incorporation was 
assessed by spectrophotometry (NanoDrop ND-1000, 
Thermo Scientific, USA) and fluorescent gel electrophor- 
esis. Three hundred nanograms of each Cy3-labelled test 
sample was competitively hybridised with 300 ng Cy5- 
labelled common reference pool on a 15 K feature custom 
microarray, following the manufacturer's instructions 
(Agilent Technologies, UK). The hybridisation reactions 
were incubated at 65°C with 10 rpm rotation for 17 hours 
in an Agilent rotary hybridisation oven and then washed 
with Gene Expression Wash Buffers 1 and 2, with a final 
wash using Stabilisation and Drying solution, again follow- 
ing the manufacturer's instructions (Agilent Technologies, 
UK). The hybridised microarrays were scanned using an 
Axon Genepix 4200A scanner with Genepix Pro 6.1 image 
acquisition software (Molecular Devices, UK) using 40% 
laser power, 5 |im pixel size resolution and auto photo- 
multiplier tube (auto-PMT) function with 0.05 saturation 
tolerances. The raw microarray images were processed 
using Agilent Feature Extraction (FE) software version 
9.5.3.1 that performed feature grid alignment, extraction 
and quantification. The fluorescence intensity results fOes 
from the FE software were imported into the GeneSpring 
GX version 12 software (Agilent Technologies, UK) for 
differential gene expression analysis. Data were normalised 
using Lowess normalisation of log2-expression ratios 
without baseline transformation. Features showing low 
quality according to Agilent quality control metrics were 
discarded to provide a final feature set for analysis. Details 
of microarray experiments 1 and 2 have been submitted to 
ArrayExpress and assigned accession numbers E-MTAB 
-1484 (Experiment 1) and E-MTAB- 1478 (Experiment 2). 
The recording of the microarray experimental metadata 



complies with Minimum Information About a Microarray 
Experiment (MIAME) guidelines. 

RT-qPCR 

To validate gene expression results from microarray 
experiments, the abundance of selected differentially 
expressed transcripts was determined by reverse transcrip- 
tion quantitative PGR (RT-qPCR). Per experiment, three 
targets that showed stable expression levels in microarray 
hybridisations were selected as reference genes. (Experi- 
ment 1: 60S ribosomal protein S20, 40S ribosomal protein 
L44 and RMD-5 homolog; experiment 2: Hypoxanthine- 
guanine phosphoribosyltransferase (HGPRT), RMD-5 
homolog and Elongation factor la). For each target 
sequence, primers were designed with a melting tempe- 
rature (Tm) of ~60°C using Primer 3 software (Additional 
file 5). Aliquots (1 [ig) of total RNA samples previously 
used in microarray analyses were reverse transcribed 
(Superscript III, Invitrogen, UK) using random hexamers 
and anchored oligo-dT in a 3:1 molar ratio. No-template 
controls and controls omitting RT enzyme were included 
on each assay plate to detect potential DNA contamin- 
ation. A cDNA pool containing equal amounts of all sam- 
ples was made and included on each assay plate, serving 
as a calibration sample (20-fold dilution) and for deriv- 
ation of a standard curve from serial dilutions. RT-qPCR 
reactions were performed in duplicate in a total volume of 
20 \iL containing 5 |iL sample cDNA (20-fold dilution), 
0.3 |iM of each primer and 10 i^L Absolute SYBR Green I 
mix (ThermoFisher Scientific, UK), using the Mastercycler 
ep realplex (Eppendorf, UK) with the following amplifica- 
tion conditions: 95°C for 15 minutes, followed by 40 cycles 
of 94°C for 30 seconds, 15 seconds at the specific primer 
pair annealing temperature (Ta; Additional file 5) and 72°C 
for 30 seconds. After amplification a melt curve from 55°C 
to 95°C at 0.5°C increments for 15 seconds each was 
performed to ensure that a single product was amplified 
in each reaction. Threshold cycles were analysed using the 
PGR cycler software. Standard curves were derived from 
plots of the threshold cycle against the logarithm of the 
relative concentration of cDNA pool. Primer efficiency (E) 
was derived from linear fits to the standard curve 
according to the equation E = j^q*"^''*'"'"^'. The BestKeeper 
tool [57] was employed to analyse expression stability of 
three reference genes and determine a robust BestKeeper 
expression index as a geometric mean for the three refer- 
ence genes, which was in turn used to establish relative 
gene expression ratios using the AACt method using the 
Relative Expression Software Tool (REST) Multiple Con- 
dition Solver (MCS) [58]. 

Statistical analysis 

Microarray gene expression data were analysed using 
GeneSpring GX version 12 (Agilent Technologies). The 
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analysis of constitutive differential gene expression in ex- 
periments 1 and 2 used Student's i-test adapted for sam- 
ples with unequal variance (Welch) using a fold change 
threshold of 1.3. The analysis of differential gene expres- 
sion induced in microarray experiment 2 employed two- 
way ANOVA to compare exposure of both salmon louse 
strains at two time points against control conditions. 
Multiple testing corrections were not applied to any stat- 
istical analysis of this gene expression study as this can 
often be over-conservative when studying potentially 
subtle gene expression responses to stimuli [56,59]. This 
decision is supported by confirmation of differential 
expression by RT-qPCR in the current study. Network 
analysis of microarray experiment 2 expression data was 
performed using the BioLayout Express^° application 
[22]. A network graph was constructed using the Pear- 
son correlation coefficient (threshold of 0.94) to deter- 
mine similarities between expression profiles, which 
were then arranged into groups of features with similar 
profiles using the Markov clustering algorithm (MCL) 
with the default inflation setting (2.2) for optimal clus- 
tering. Gene enrichment analysis was performed on lists 
of features chosen based on differential gene expression 
patterns using default settings of the FuncAssociate 2.0 
web application [60]. Gene enrichment was calculated 
according to the significance (p< 0.05) of the association 
between the list of features and the GO attributes repre- 
sented on the microarray. Relative expression ratios 
from RT-qPCR experiments were tested for normality 
and equal variance and log transformed to allow as- 
sumptions to be satisfied before being subjected to one 
way ANOVA using Minitab 16.1 software (Minitab Inc., 
UK). The significance level was set at p<0.05 in all tests. 
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